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Abstract As the largest mass concentrations in the local 
Universe, nearby clusters of galaxies and their central galax- 
ies are prime targets in searching for indirect signatures of 
dark matter annihilation (DMA). We seek to constrain the 
dark matter annihilation emission component from multi- 
frequency observations of the central galaxy of the Virgo 
cluster. The annihilation emission component is modeled 
by the prompt and inverse-Compton gamma rays from the 
hadronization of annihilation products from generic weakly 
interacting dark matter particles. This component is fitted 
to the excess of the observed data above the spectral en- 
ergy distribution (SED) of the jet in M87, described with 
a best-fit synchrotron-self-Compton (SSC) spectrum. While 
this result is not sufficiently significant to claim a detection, 
we emphasize that a dark matter "double hump signature" 
can be used to unambiguously discriminate the dark matter 
emission component from the variable jet-related emission 
of M87 in future, more extended observation campaigns. 

Keywords dark matter ■ clusters of galaxies ■ high-energy 
emission ■ extragalactic jets ■ M87 



^ 1 Introduction 
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The thermal freeze-out of weakly interacting particles (so- 
called 'WIMPs') with masses at the electroweak symme- 
try breaking scale ZSgw = 1/ V2'^^Gf ~ 246 GeV leads 
to a relic density agreeing with the observed one Qdmh^ - 
0.1123 + 0.0035 |1,2|. This so-called WIMP miracle pro- 
vides a natural solution to the dark matter problem. In the 
present-day Universe, WIMPs can still annihilate in regions 
with large mass concentrations. The Galactic Center iT"?!, 
dwarf galaxies in the halo of the Milky Way tSj^J^^^, and 
galaxy clusters HlOlll II have been studied to obtain Umits on 
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dark matter annihilation emission. Annihilation gives rise to 
pairs of heavy quarks, leptons or vector bosons, which decay 
leading to the emission of gamma rays, electrons, positrons, 
and neutrinos. The secondary radiation shows a peak at £'ew/20 
10 GeV reflecting the high multiplicity of secondary par- 
ticles in a typical annihilation event. The canonical peak 
energy lies at the high end of the bandwidth of the Fermi- 
LAT detector, and at the low end of the accessible energy 
range for ground-based Cherenkov telescopes. A very well 
motivated WIMP candidate is the Lightest Supersymmet- 
ric Particle (LSP): the neutralino fTZ]. Annihilating neutrali- 
nos produce vector bosons, leptons or quarks. The differen- 
tial gamma-ray flux |fT3]| produced through their subsequent 
hadronization and decay from a source at distance D and 
volume V is given by: 
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where / is the boost factor that accounts for enhancement 
due to sub-halo clumping of dark matter in the M87 halo, 
(cr^v) is the thermally averaged annihilation cross section, 

is the mass of the dark matter particle, and dNy/dE is the 
gamma photon spectrum coming from the decay of neutral 
pions from hadronization in the annihilation process (prompt 
pion emission) lfT4ll . p is taken to be pnfw - PNPw(r) the 
Navarro-Frenk- White dark matter density profile fl5\ of the 
dark matter halo obtained from numerical simulations. 

A minimum mass scale of 10"^ for the sub-halo clum- 
ping has been inferred from the transfer function of density 
perturbations in the early Universe for weakly interacting 
dark matter particles |16|. However, tidal interactions with 
the baryon-dominated cores of dark matter halos and super- 
nova feedback could destroy these structures to a large ex- 
tent. Therefore, the boost factor / = / p^dV/ J p^p^dV is 
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introduced as a free parameter to account for the unknown 
enhancement due to sub-halo clumping. 

The inverse-Compton interaction is the up-scattering of 
photons by high-energy charged particles. The differential 
gamma-ray flux given by Eq. (2) describes inverse-Compton 
scattering off' the cosmic microwave background by rela- 
tivistic electrons and positrons. b(E') is the total rate of 
electron/positron energy loss due to inverse-Compton scat- 
tering as in |i,17J . P{E,E') is the differential power emitted 
into photons of energy E by an electron/positron with en- 
ergy E', and dNe/dE is the spectrum of secondary electrons 
and positrons fT4\ : 
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The main challenge in constraining the putative dark mat- 
ter annihilation component is to discriminate it against gam- 
ma rays from astrophysical sources and cosmic ray interac- 
tions. Here, we show as an exemplary case study how recent 
multi-frequency data of the center of the Virgo cluster, har- 
boring the giant cD galaxy M87 with its gamma-ray emitting 
radio jet, can be used to constrain dark matter particles. Note 
that we do not refer to the inverse-Compton scattering in so- 
called leptophilic dark matter annihilation scenarios, but in- 
clude the inverse-Compton emission component due to the 
electrons and positrons from the decay of isospin-symmetric 
annihilation products. 

We adopt a distance of (16.5 + 1) Mpc [!l8| to M87 and 
therefore omit redshift corrections to the energy through- 
out the paper Section 2 describes the data sets chosen for 
the study. The radiation code employed to model the spec- 
tral energy distribution of M87 is explained in Section 3, 
and physical parameters inferred from the fit of the data are 
briefly discussed to show that the model is a viable interpre- 
tation of the nonthermal particle content in the jet of M87. 
The generic model for the weakly interacting massive parti- 
cles and their radiative signatures is described in Section 4, 
and a fit of a component from dark matter annihilation to the 
excess above the astrophysical model is presented. Finally, 
we discuss the results and draw conclusions regarding future 
search strategies. 



2 Data 

The dataset used here is comprised of observations of M87 
by the Chandra X-ray Observatory, Fermi-LAT, MAGIC, 
and the H.E.S.S. system of Cherenkov telescopes 11911201 
I2ni22ll . As presented in |20|, additional observational data 
exist for the radio-to-optical regime. However, we do not 



use these to constrain the fit, but rather only require the pro- 
jected emission to not exceed these observations, since in 
this wavelength regime the unavoidable contamination due 
to starlight and dust from the central region of M87 is very 
difficult to assess. There may also be hidden nonthermal 
components in this energy region unrelated to the emission 
region where the high-energy emission originates from. 

The data sample was not taken contemporaneously, but 
we carefully checked to avoid inclusion of an observation 
containing a significant flare. All the data used here can thus 
be considered a representative long-term average, reason- 
ably well describing the steady and low-state spectral energy 
distribution of the high-energy emission component in M87. 



3 Synchrotron-self-Compton emission 

The nonthermal emission from the relativistic jet emerging 
from the nucleus of M87 follows a spectral energy distribu- 
tion across almost 20 orders of magnitude which can be de- 
scribed by the synchrotron and self-Compton radiation pro- 
duced by electrons (and positrons) accelerated in the jet, pre- 
sumably by shock waves. Whereas the high amplitude flares 
track single shock waves or magnetic reconnection events 
in the expanding flow, guided by the helical structure of the 
magnetic field, the steady-state (long-term average) emis- 
sion which dominates the energy output from the jet of M87 
corresponds to the superposition of a large number of shock 
waves 1231 . a strong stationary shock such as the reconfine- 
ment shock Il24l or the emission from a sheath surrounding 
the spine of the jet ll25l . 

To obtain a fit for the broad-band spectral energy dis- 
tribution of M87, we used an implementation of the SSC 
model 1261 in which the cooUng break in the electron spec- 
trum is self-consistently determined, and the cross section 
for Compton scattering in the Klein-Nishina regime is ac- 
curately treated. We also took into account the inhomoge- 
neous nature of the emission zone, by considering only the 
emission from the inner jet for the fit, i.e. we considered the 
radio to optical emission, related to the emission from larger 
scales in the jet, as an upper bound for the model. 

Adopting for the Doppler factor 6-3.9 (bulk Lorentz 
factor r - 2.3) as in the model fitted to the data by |20] 
we obtain a fit withx^/v = 2.5. Physical parameters of the 
emission zone are the source radius Rb = 3.5 x lO'^* cm 
and the magnetic field strength B = 3 G. For the injected 
electron power law distribution with exponential cut-off, we 
obtain a differential slope of s = 2.2, maximum Lorentz fsic- 
toT r„ax - 10^ and normalization factor = lO^cm^^s^^ 
(cf. Il26l ). This results in an injection luminosity of L,„j = 
3 X 10^' erg s"' consistent with the energetics of the jet in- 
ferred from its large scale radio structure |27|. Model fits 
including the low-energy continuum l,20J yield different re- 
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Fig. 1 Result of the combined fit including the spectral energy distri- 
butions due to the SSC mechanism in the jet of M87 and the annihila- 
tion of dark matter particles. The boost factor / accounts for sub-halo 
clumping, and m denotes the generic WIMP mass. 

suits, but fail to produce an acceptable fit including the very 
high energy observations. 

For the black hole mass of M87 Mbh = 6.4 x IO^Mq, 
the Eddington luminosity is given by Le = 8.32 x 10^' 
erg s"' . Thus, the nonthermal power release amounts to only 
~ IO^^Le- Although a sub-Eddington state of the black hole 
is generally expected for high-peaked blazars, the extremely 
low power is peculiar In fact, it is not possible to model the 
data by simply adjusting the Doppler factor of an SSC fit 
obtained for high-peaked blazars, i.e. by changing the incli- 
nation of the jet axis with respect to the observer The idea 
of treating M87 as a misaligned blazar Il28l does not seem 
to be sufficient. Furthermore the large difference between 
break energy and cut-off is not typical for blazar spectra. 
This might be related to the extremely low accretion rate, 
indicating that the AGN is fading out due to a lack of acc- 
retable matter, or to the fact that the observed emission is 
dominated by the emission from a low-bulk-Lorentz-factor 
sheath surrounding the jet. 

The SSC fit of the data obtained in this way provides an 
accurate and practically unique model for the SED, and the 
small size of the emission region is in line with observations 
of short-time variability. 

4 Fitting the excess emission component with dark 
matter annihilation 

To now fit an additional component due to dark matter an- 
nihilation, and study whether the statistical agreement be- 
tween model and observations can thereby be further im- 



proved, we assume a generic species of annihilating WIMPs 
with rest mass in the GeV-TeV range and a thermally aver- 
aged annihilation cross section of (cr^v) = 3 x 10"^"^ cm^s"'. 
It has been shown that many models with multi-TeV masses 
provide large pair annihilation cross sections which are still 
in agreement with the thermal freeze-out of these particles in 
the early Universe ||29]| . The choice of the annihilation cross 
section also reflects the possible existence of boost factors 
in the particle physics sector and should be considered as 
moderate upper limit (cf ll30l ). For the dark matter distribu- 
tion of M87, we use the analysis of [31], resulting in a de- 
scription of the halo according to the Navarro-Frenk- White 
flS I model. The normalization of the intensity of the result- 
ing emission is furthermore fixed by choosing the boost fac- 
tor from unresolved substructure in the halo of M87, which 
from recent numerical experiments |32J is expected to be of 
order / = 10^-10^. The spectra for the emission due to the 
decay of charged and neutral pions from hadronization in 
the annihilation process are generated using the DarkSUSY 
code lf33l . A total of 10^ realizations of neutrahnos are pro- 
duced to obtain the average secondary spectra for generic 
WIMPs motivated by supersymmetric theories. 

Here we also include a treatment of the inevitable con- 
tribution of inverse-Compton emission from energetic elec- 
trons/positrons from the decay of the charged pions up-scat- 
tering cosmic microwave background photons. By using the 
known number density of these 2.7 K background photons 
(413 cm"-'), the respective contribution from inverse-Comp- 
ton emission can be of the same order as the prompt pion de- 
cay emission. This results in a telltale "double hump struc- 
ture" of the dark matter related emission. Employing a x^- 
test, we search for the best-fit model of the active galactic 
nucleus (AGN) as discussed in the previous section, and the 
complete dark matter related model. The results of the fit are 
shown in Fig. 1. 

A minimal value of 1 .6 is achieved for a particle mass 
of 4.7 TeV and a boost factor of 812, both well within the 
range of values accessible in state-of-the-art particle physics 
and numerical models, as discussed in the previous section. 
Assuming that the main multi-wavelength features of such 
a generic dark matter model can be described by two addi- 
tional parameters (neutralino mass and boost factor) com- 
pared to the SSC model alone, it can be shown by nonlinear 
regression analysis that the deviation of the data points with 
respect to the SSC model amounts to 2cr. We note that this 
result is not significant enough to claim the realization of a 
certain dark matter model, but even our simplified assump- 
tions show the capability of multi-wavelength methods for 
unvealing the characteristics of the constituents of dark mat- 
ter In Fig. 2, we show the resulting spectral energy distribu- 
tion from the combined SSC-DMA model. 
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Fig. 2 Spectral energy distribution of the best-tit combined SSC model for M87 and dark matter model, including all data points used in this 
analysis 1 1912012 11221 . The inset legend provides the mass and boost factor for the generic WIMP with cross section {(Tav) = 3 x 10"^'* cm^s"' . 



5 Discussion and conclusions 

We show that while synchrotron-self-Compton emission can 
fit the observations reasonably well, the introduction of a 
component due to annihilating dark matter particles margi- 
nally improves the agreement between theoretical modeling 
and multi-frequency observations of the long-term average 
SED of the radio galaxy M87. It is essential to include the 
inverse-Compton emission component in the dark matter an- 
nihilation model. The method thus shows great sensitivity 
in ruling out dark matter models at the TeV scale which is 
innate to some supersymmetric extensions of the Standard 
Model, while annihilation cross section and additional boost 
from substructures are concordant with the paradigm of ther- 
mal freeze-out of such a particle at the onset of hierarchical 
structure formation. Although the significance of this result 
is not sufficient to claim evidence of a specific particle, it is 
encouraging further studies. 

Since the SSC component generally shows variabihty 
while the DMA component remains steady, we emphasize 
that the significance of the marginal excess reported here 
would increase during low states of the SSC emission. There- 
fore we encourage to extend the temporal coverage of multi- 
wavelength observations of M87. Measuring the putative 
dark matter "double hump structure" in the spectrum with 



a high significance would open up the possibility to extract 
the dark matter properties (WIMP mass, annihilation cross 
section and boosting due to substructure) more accurately. 

The high mass scale of the generic WIMP favored in 
our model here is in-line with recent findings at the LHC 
||34l . A high mass scale may also help in explaining the 
measured extragalactic gamma-ray background (EGB) 135)1 , 
after a careful reassessment of the Fermi-measured EGB 
ll36l . Detection of high-mass WIMPs by elastic scattering 
experiments [371 will be difficult due to the suppression of 
the event rate by the low number density of WIMPs in the 
Galactic halo. 
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